subroutine var_init
use storage
implicit none
integer i_x,i_y,i_z
real(myk) :: temp
real(myk), dimension(1-ng_x:nt_x+ng_x) :: temp_rho,temp_rhu,temp_rhe,temp_rya,temp_temp
open(unit=100,file='1d_flame.dat')
   read(100,*)
   read(100,*)
   do i_x=1,nt_x
      if((i_x.ge.istart_x-ng_x) .and. (i_x .le. iend_x+ng_x)) then
         read(100,*) temp,temp_rho(i_x),temp_rhu(i_x),temp_rhe(i_x),temp_rya(i_x),temp_temp(i_x)
      else 
         read(100,*) temp,temp,temp,temp,temp,temp
      endif
      if(i_x == 1) then
         temp_rho(1-ng_x:0) = temp_rho(i_x)
         temp_rhu(1-ng_x:0) = temp_rhu(i_x)
         temp_rya(1-ng_x:0) = temp_rya(i_x)
         temp_rhe(1-ng_x:0) = temp_rhe(i_x)
      endif
      if(i_x == nt_x) then
         temp_rho(nt_x+1:nt_x+ng_x) = temp_rho(i_x)
         temp_rhu(nt_x+1:nt_x+ng_x) = temp_rhu(i_x)
         temp_rya(nt_x+1:nt_x+ng_x) = temp_rya(i_x)
         temp_rhe(nt_x+1:nt_x+ng_x) = temp_rhe(i_x)
      endif
   enddo
close(100)
print*,'1d check rho: ',temp_rho
print*,'1d check rhu: ',temp_rhu
print*,'1d check rhe: ',temp_rhe
print*,'1d check rya: ',temp_rya
!read(*,*)
do i_z=istart_z-ng_z,iend_z+ng_z
   do i_y=istart_y-ng_y,iend_y+ng_y
      do i_x=istart_x-ng_x,iend_x+ng_x
         var_rho(i_x,i_y,i_z) = temp_rho(i_x)
         var_rhu(i_x,i_y,i_z) = temp_rhu(i_x)
         var_rhe(i_x,i_y,i_z) = temp_rhe(i_x)
         var_rya(i_x,i_y,i_z) = temp_rya(i_x)
      enddo
   enddo
enddo
do i_z=istart_z-ng_z,iend_z+ng_z
do i_y=istart_y-ng_y,iend_y+ng_y
do i_x=istart_x-ng_x,iend_x+ng_x
   var_rhv(i_x,i_y,i_z) = 0._myk
   var_rhw(i_x,i_y,i_z) = 0._myk
enddo
enddo
enddo
end subroutine var_init
